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A Dzyaloshinskii-Moriya Anisotropy in Nanomagnets with in-plane Magnetization 
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We report on a new source of in-plane anisotropy in nanomagnets due to the presence of 
Dzyaloshinskii-Moriya interaction (DMI). This anisotropy depends on the shape of the magnet, 
and is orthogonal to the demagnetization shape anisotropy. This effect originates from the DMI 
energy reduction due to an out-of-plane tilt of the spins at edges oriented perpendicular to the mag¬ 
netization. Our investigation combining experimental, numerical and analytical results demonstrate 
that this energy reduction can compensate the demagnetization shape anisotropy energy in magnets 
of elongated shape, provided that their volumes is small enough and thus that their magnetization 
is quasi-uniform. 


The Dzyaloshinskii-Moriya interaction (DMI), which 
is the antisymmetric analog of the Heisenberg exchange 
energy, reduces the energy of a magnetic system when 
neighboring spins Si and Sj are not parallel, and can be 

expressed in an Hamiltonian as dij • (^Si x Sj^, where 

dij is the Dzyaloshinskii-Moriya vector. The DMI is 
known to be a direct manifestation of spin-orbit coupling 
(SOC) in systems with broken inversion symmetry [1, 2]. 
During the last couple of years, DMI has triggered an 
increasing interest and a lot of research in magnetism 
community because it could lead to several technological 
breakthroughs in the area of magnetic recording. In par¬ 
ticular, the DMI can impose the structure of magnetic 
domain walls (DW) in perpendicular magnetic films to 
change from conventional Bloch into chiral Neel configu¬ 
rations [3-5], controlling the direction of propagation and 
increasing dramatically the velocity of DW in racetrack 
structures [6-8]. Interestingly, this chiral interaction can 
also lead to new magnetic textures, such as spin spirals 
or magnetic skyrmions of a few atomic lattice periodic¬ 
ity/size [9, 10] that hold promises as a new approach for 
magnetic memories [II-I3]. 

Large interfacial DMI of magnitude comparable with 
the Heisenberg exchange can be brought about in thin 
metallic multilayers made of 3d ferromagnetic {e.g. Fe, 
Co, Ni) and bd heavy metals {e.g. Ta, W, Ir, Pt), in 
which the inversion symmetry is broken at interfaces 
and the large SOC is provided by the heavy metal in 
contact with the ferromagnetic layer [6, 7, 9, 12, 14]. 
The impact of the DMI in magnetic materials with out- 
of-plane magnetization were experimentally observed by 
several means, such as the measurement of the asymmet¬ 
ric growth and propagation of magnetic domains or DW 
with external magnetic fields [6, 15-18], the tilting of the 
magnetic DW in the presence of an external transverse 
magnetic field [19], the local measurements of the mag¬ 
netization or its stray field using scanning probes [5, 20], 
the study of the spin wave propagation in magnetic mul¬ 
tilayers with Brillouin light spectroscopy [21], or the ob¬ 
servation of magnetic skyrmions by spin-polarized STEM 


or scanning transmission x-ray microcopy with polarized 
light [22-24]. Most of these techniques and observations 
rely on the particular structure of the DW induced by 
the DMI, namely a Neel configuration in which the chi¬ 
rality (sense of rotation of magnetization along an axis 
lying in the DW plane) is fixed by the sign of the DMI. A 
DMI-induced tilt of the magnetization close to the edge of 
out-of-plane magnetic structures has also been predicted 
[13, 25], but there are no experimental evidences of this 
phenomenon yet. 

For symmetry reasons, most of the studies of the DMI 
have been performed with perpendicularly magnetized 
ultra-thin films. In this letter, we investigate the effect 
of the DMI in magnetic films with in-plane magnetiza¬ 
tion. We demonstrate the existence of an original in¬ 
plane anisotropy term due to the tilt of the magnetiza¬ 
tion that is induced by the DMI at edges of magnetic 
nanostructures [25]. We first present experimental ob¬ 
servations of the unexpected magnetization direction in 
sub-micrometer-scale ellipses made of a thin layer of in¬ 
plane magnetized CoFeB sandwiched between a Pt and 
an MgO layer: For the smallest observed structures, the 
magnetization lies along the minor axis of the ellipses. 
This behavior is well reproduced in micromagnetic sim¬ 
ulations in which the interfacial DMI is considered. This 
observation also allows the micromagnetic DMI ampli¬ 
tude D to be indirectly determined, providing that other 
micromagnetic parameters are known. We then explain 
the origin of the effect and give an approximated analyt¬ 
ical formulation for the critical D value, D±, that brings 
the magnetization from lying along the major axis to the 
minor axis in rectangular prisms. Finally, we present a 
set of phase diagrams of the magnetic configurations as 
a function of these micromagnetic parameters. 

The studied samples are composed of 
MgO|Co 2 oFe 6 oB 2 o (1-5 nm)|Pt (0.8 nm) stacks that 
are grown by sputtering deposition. We fabricated elon¬ 
gated nanostructures (nominally ellipses) with variable 
areas and aspect ratios, using e-beam lithography and Ar 
ion beam etching. The major axis of the studied ellipses 
nominally ranges from I jam to 100 nm, and their aspect 
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Figure 1. (Color online) AFM/MFM measurement and ge¬ 
ometry of the problem, (a) Actual measurement showing in 
pseudo-3D the topographical signal (AFM) colored with the 
phase of the MFM signal. The MFM phase is coded blue 
or red for attractive or repulsive forces, allowing the mean 
magnetization to be determined, as indicated by the arrow. 
(b,c) Geometries used to perform micromagnetic numerical 
simulations and the analytical calculations (c). 

ratios from 1 to 10 [Fig. l(b,c)]. We have used atomic 
and magnetic force microscopy (AFM/MFM) to image 
both the topographies (real sizes and shapes) and the 
magnetic configurations (direction of the magnetization) 
of these nanostructures. Particular attention should be 
taken about the choice of a tip with ultra-low moment 
to avoid perturbing the sample magnetization by the 
tips stray field. A typical elliptical structure is displayed 
in Fig. 1(a) in which the topographical signal (AFM) 
is rendered in pseudo-3D and the MFM phase signal is 
color-coded. 

The experimental measurement of a series of ellipses of 
different sizes is presented in Fig.2(a). It illustrates how 
the magnetization turns from the major axis for large 
ellipses with L « 0.7 jam (left MFM image) to the mi¬ 
nor axis for a small ellipse with L ^ 160 nm (right MFM 
image). In the intermediary case we found one ellipse 
{L ^ 250 nm) with the magnetization pointing neither 
along the minor nor the long axis. As we will see below, 
this state is probably a limit case where the DMI-induced 
anisotropy equals the shape anisotropy. We emphasize 
that such effects were not observed on a symmetric ref¬ 
erence sample made of MgO|Co 2 oFe 6 oB 2 o (2.9nm)|MgO. 

We reproduce this behavior using micromagnetic simu¬ 
lations, for which we used the GPU-based MuMax3 code 
(version 3.6.1 and 3.8) which natively includes the in¬ 
terfacial DMI [26]. In our simulations, as the simulated 
systems are thinner than the exchange length, we con¬ 
sider that the magnetization is uniform over the thick¬ 
ness. The system is hence modeled with a single layer, 
and consequently discretization cells have their size along 
z determined by the thickness t, while sizes of 0.4 to 
0.8 nm were considered along x and y directions. In order 
to define the geometrical shape of the simulated struc¬ 
ture, we have used the exact sample dimensions deter¬ 
mined by the AFM characterization (taking in account 
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Figure 2. (Color online) DMI-induced anisotropy orthogonal 
to the shape anisotropy as observed by MFM and simula¬ 
tions. All scale bars are 200 nm long. The arrows indicate the 
mean magnetization, (a) Experimental MFM phase images 
are displayed. The actual topographical shape is indicated 
by a dotted line. The MFM phase is coded in the same way 
than in Fig. I. (b) Corresponding simulations for different D 
values (all the other micromagnetic parameters fixed), using 
an ideal dipolar MFM tip 30 nm above the magnetic layer. 
The observations of the magnetization direction in ellipses of 
different shapes (but same materials) allow constraining the 
Dzyaloshinskii-Moriya interaction strength D: The shadings 
indicate the D values compatible with observations. In the 
intermediate case, the magnetization does not align along one 
of the axes, reflecting the minute energy difference between 
the states with different magnetization directions, (c) Map of 
the DMI and exchange energy densities (darker means lower 
energy density) for the case D = 2.1 mJ/m^. The energy gain 
concentrates close to the edges. 


the effect of the tip width). The magnetic parameters 
used for the simulations are the magnetization at satura¬ 
tion Ms = 1.06 ± 0.08 MA/m and the crystalline per¬ 
pendicular anisotropy Ku = 0.26 ± 0.06 MJ/m^ both 
determined by SQUID measurements [27], and an ex¬ 
change stiffness A of 10 pJ/m typically found in the liter¬ 
ature for this material [28]. The value of the micromag¬ 
netic DMI parameter D is still not available for a large 
range of interfaces, but a rather large value is expected 
at the interface between Pt and ferromagnetic materials 
[4, 6, 7, 16, 18, 21]. Simulations were hence performed 
for different values of D and, as displayed in Fig. 2(b), 
there is a change of magnetization direction as the DMI 
increases. Comparing simulations with MFM measure¬ 
ments allows the estimation of D ^ 2.llQ;3mJ/m^ at 
the CoFeBjPt interface. This value can be compared. 
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for example, to the Brillouin light scattering experiments 
where D = 1.2mJ/m^ for the Co(1.2nm)|Pt interface 
[21] or 0.8 for CoFeB(0.8 nm)|Pt [29], or to domain wall 
experiments in Co(0.6 nm)|Pt where I) = 2.2mJ/m^ [18]. 
It is important to note that the value of D is very sen¬ 
sitive to the actual magnetic parameters of the material 
(as readily visible on the phase diagrams of Fig. 3, e.g. 
D ^ 2.4mJ/m^ if A = 15pJ/m is used). However, pro¬ 
vided a precise determination of the micromagnetic pa¬ 
rameters (Mg, Ku and H), the observation of nanoscale 
ellipses could hence lead to a novel and more direct way 
to estimate D. One should finally note that the energy 
minimization of the simulations must be studied with 
caution because the energy landscape is extremely flat 
close to the transition. This also explains why, in the 
real system, minute inhomogeneity might stabilize the 
magnetization along any direction (see for example the 
intermediate case of Fig. 2). 

In order to get a better insight into the mechanisms 
at play, an analytical model can be formulated. In mag¬ 
netic thin films with moderate DMI, the DMI induces a 
tilt of the magnetization at the border. In perpendicu¬ 
larly magnetized films, this tilt occurs along the whole 
border of the structure [25], while for in-plane magne¬ 
tized films with interfacial induced DMI, the tilt only 
occurs at the borders with an orthogonal component of 
the magnetization. In particular, it is important to no¬ 
tice that there is no tilt if the magnetization is parallel 
to the edge, while it is maximal for a magnetization di¬ 
rection perpendicular to the border. For the latter case, 
the tilt of magnetization reduces the system energy by 
[25, 30]: 


^tilt — — 






--tL 


where L is the border length, t is the thickness of the 
magnetic film. The parameter Ky^ is the effective pla¬ 
nar anisotropy constant and is expressed as Kyz = 
(A4 —My) — Ku for a simple shape with a remain¬ 
ing crystalline perpendicular anisotropy Ku^ and Mi are 
the demagnetizing factors. As the borders where tilting 
occurs depend upon the direction of the in-plane mag¬ 
netization, this phenomenon induces a planar anisotropy 
that competes with the dipolar-induced shape anisotropy. 

To describe this effect, let us consider a simple shape: a 
rectangular prism with dimensions L, W and t [along the 
X, y, and z axes respectively as shown in Fig. 1(c)]. We 
consider a system small enough to be in a single magnetic 
domain state. The energy difference Ex — Ey between the 
states with the magnetization along x and y is then given 


by: 


Ex ~ Ey 
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where Kxy = (A4 — My) is the in-plane shape 

anisotropy and Kxz = (A4 — Mx) — Ku- This 

DMI-induced tilt results in an in-plane anisotropy of di¬ 
rection orthogonal to the shape anisotropy, i.e. favor¬ 
ing the magnetization parallel to the smaller side of the 
structure. In the case of a rectangular parallelepiped, the 
“DMI-anisotropy” can be approximated as: 


A'dmi 


L-W 

2^AKMz LW 


The threshold value at which the magnetization 
switches from one axis to the other is then: 


Dl «(A4 - , (1) 


which is composed of a part that depends on shape and 
volume and another part that depends almost only on 
the material magnetic parameters. 

To illustrate this effect more comprehensively, in Fig. 3, 
using micromagnetic simulations, we map out phase dia¬ 
grams for elliptical shapes [see Fig. 1(b)] of the magneti¬ 
zation direction as a function of the magnetic parameters 
Ku^ Mg, A and the geometrical parameter t using typ¬ 
ical parameters found for thin magnetic films. In each 
phase diagram, three regimes that are controlled by the 
DMI are observed: (I) at low DMI, the magnetization lies 
along the major axis as expected in classical systems, (2) 
at larger values of DMI, the DMI-anisotropy overcomes 
the shape anisotropy and aligns the magnetization along 
the minor axis, and (3) for very large DMI, the magneti¬ 
zation is not anymore single domain, but multi-domains, 
helices or skyrmions appear, because the presence of do¬ 
main walls minimizes the energy, as each of these DW, 
being of Neel type, reduces the DMI energy. Illustra¬ 
tions of typical configurations are given in the top part 
of Fig. 3. The tilt at the edges perpendicular to the mag¬ 
netization are also clearly visible (color-coded black and 
white for down or up magnetization). 

The analytical model (Eq. I) can be used to plot D_\_ 
in a simplified case where the demagnetization factors of 
the rectangular prism are considered for a uniform mag¬ 
netization [31]. As shown in Fig. 3, the analytical curves 
are reproducing the trends of the micromagnetic simu¬ 
lated system. The difference between D± in the micro- 
magnetic simulations (the border between the red right 
triangles with the green up triangles) and the analytic 
curve cannot be simply explained by the Mi of the differ¬ 
ent shapes {i.e. demagnetizing factors for rectangles and 
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Figure 3. (Color online) Numerical and analytical phase di¬ 
agrams. Simulated state (ellipses) of for different Ku (a), 
Ms (b), A (c) and t (d) as a function of D, and the analyti¬ 
cal critical D± for rectangles. Right pointing triangles (red) 
correspond to states with the magnetization along x, up tri¬ 
angles (green) along y, spheres (blue) are used if the state 
is very inhomogeneous with several domains, and hexagons 
(open) if the state is with perpendicular magnetization. The 
physical parameters of the simulations are Ms = IMA/m, 
Ku = 377kJ/m^ A = 15pJ/m, L = 210, VF = 70 and 
t = 1.5 nm if not indicated otherwise. Four typical micro- 
magnetic states of panel (a) which correspond to the black 
symbols are drawn on top of the figure. 
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Figure 4. (Color online) Analytical critical D for trans¬ 
verse magnetization as a function of the long size of a rect¬ 
angle for different aspect ratio and Ku- The uniaxial perpen¬ 
dicular anisotropy plays a crucial role and diminishes greatly 


ellipses). Indeed, simulations were also realized for rect¬ 
angular nanomagnets showing a behavior much closer to 
the ellipses than to the analytical model. The compar¬ 
ison of the different components of the energy indicates 
that the main limitation of the analytical model (which 
uses a strictly uniform magnetization) is that the de¬ 
magnetization energy is very much reduced by the tilting 
occurring at the edges, hence reducing the needed DMI- 
anisotropy necessary to align the magnetization along the 
minor axis/short side of the structures. Nevertheless, the 
simple analytical model allows a very quick assessing of 
the potential role of the DMI-anisotropy in a system. 

In Fig. 4, the dependence of the analytical critical D±^ 
for transverse magnetization is plotted as a function of 
the length (L) for different aspect ratios (L/W) and Ku- 
A larger value of the out-of-plane anisotropy Ku favors 
the DMI-anisotropy, because the tilting at edges increases 
(in absolute value) the associated negative energy term. 
As expected, D± decreases with increasing aspect ratio, 
and smaller structures require weaker DMI to induce the 
transverse stabilization of the magnetization. In partic¬ 
ular, it is interesting to note that nanostructures of the 
size of typical MRAM cells will undoubtedly be sensitive 
to this effect if Pt is in contact to the magnetic layer, as 
a DMI as high as 0.5mJ/m^ are certainly expected. 

Finally, we note that the energy difference between 
states with magnetization along x or y reaches values 
greater than 1 eV, when far enough from the D± line 
on the phase diagram {e.g. \D — D±\ > 0.2mJ/m^ in 
the case of the “standard” nanomagnet of Fig. 3). Such 
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energy difference indicates a good stability against tem¬ 
perature of the states induced by the DMI-anisotropy, as 
- Ey\ /{ksT) ^ 40 for leV at 300K. 

In conclusion, we observe and explain an original 
“DMI-anisotropy” for in-plane magnetized nanomagnets. 
We show that, for small volumes and reasonable aspect 
ratios, a strong DMI-associated energy is taking over 
the shape anisotropy demagnetization energy, and that 
the magnetization is aligning along the minor axis/short 
side of elongated structures. We can note that the en¬ 
ergy difference between these states is of the order of 
the electron-volt, meaning that they are thermally very 
stable. We also present MFM observation constituting 
the first experimental evidence of the tilting of the mag¬ 
netization at the edges of magnetic structures due to 
DMI. This phenomenon used in giant magnetoresistance 
or tunnel magnetoresistance structures could possibly be 
used to electrically detect small magnetic field with a 
device of extremely small size that maybe useful, for ex¬ 
ample, in detectors or magnetic read heads. 
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